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The intrinsic optical nonlinearities of linear structures, 
including conjugated chain polymers and nanowires, 
are shown to be dramatically enhanced by the judi¬ 
cious placement of a charge diverting path sufficiently 
short to create a large phase disruption in the dominant 
eigenfunctions along the main path of probability cur¬ 
rent. Phase disruption is proposed as a new general 
principle for the design of molecules, nanowires and 
any quasi-lD quantum system with large intrinsic re¬ 
sponse and does not require charge donor-acceptors at 
the ends. © 2015 Optical Society of America 
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The design and realization of ultrafast nonlinear optical ma¬ 
terials with large responses to optical frequency electric fields 
remains an active field of pure and applied research[l-4]. To 
date, no general design rules for obtaining large nonlinearities 
from any structure have been articulated, and the response of 
materials remains well below that allowed by quantum physics. 
In this letter, we propose a general principle that may explain 
why modern molecules fall short of their potential and use it 
to demonstrate simple structures with nonlinear responses ap¬ 
proaching the physical limits. 

The off-resonance electronic nonlinear optical response of a 
molecule is completely determined by its energy spectrum and 
transition moments. Normalized to its maximum value, the 
sum over states expression for the intrinsic first hyperpolariz¬ 
ability along the x-axis may be written as [5] 
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where the sum is over all states, is the n, i element of the po¬ 
sition matrix with x = x — xqq, — Eq is the difference 

in energy of eigenstates n and 0, N is the number of electrons 


and m their mass. We include as many as 100 energy eigen¬ 
states to calculate but systems with large ^ are always well 
approximated using only three states. For brevity, we focus on 
the first hyperpolarizability, though our results also hold for the 
second hyperpolarizability, jxxxx- For the remainder of this let¬ 
ter, all hyperpolarizabilities are divided by their maximum val¬ 
ues and represent intrinsic tensor properties. It is evident from 
Eqn. (1) that a system optimized for nonlinear optics will nec¬ 
essarily have eigenfunctions with a large degree of overlap as 
well as a large change of dipole moment between contributing 
levels. This is an essential trade-off to achieve a large response. 

It is generally recognized that the hyperpolarizabilities of 
most molecules fall at least 30 times short of the limits. [5, 6]. 
Monte Carlo simulations discovered the optimum energy spec¬ 
trum for such molecules and provided a strong indicator of 
the origin of the gap [7]. Optimization of the effective poten¬ 
tial energy an electron experiences across the main direction of 
a quasi-linear molecule showed that the hyperpolarizabilities 
may be increased by tuning a few parameters [8], by modulation 
of conjugation along a chain[9], or by donor-acceptor substitu¬ 
tion and insertion of spacers[10]. However, there is no general 
rule about how to construct the ideal potential energy profile 
across a molecule to maximize the nonlinear optical response. 

In contrast, quantum graph models can sample the vast 
space of geometries and topologies of quasi-one dimensional 
structures, leading to structures whose ^ and 7 approach the 
fundamental limits [11]. A quantum graph is a collection of one¬ 
dimensional line segments, along which electrons flow freely, 
connected together in a network to constrain particle dynam¬ 
ics to one-dimensional segments within a higher dimensional 
space. Quantum graphs are well suited for exploring nonlinear 
optics because they capture the geometric and topological prop¬ 
erties of charge transfer while being exactly solvable systems 
with quasi-quadratic energy spectra. Our approach is to iden¬ 
tify the essential features of the transition moments and how 
to tailor them through topological and geometrical properties 
without the need of a potential on the graph. 

We show that quasi-linear structures, such as polymer 
chains, generate large responses when the component of the 
ground and first few excited state wavefunctions along the 
main chain direction have an abrupt phase disruption caused by 
the presence of an alternate pathway for electron flux, such as 
a short side group, or a local defect such as a Dirac delta po¬ 
tential. This effect alone generates nonlinearities near the fun- 
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Fig.l. The one-prong quantum graph. All calculations use 
a = 0.63 and c = 0.37. The results are scale-independent. 
The horizontal direction is the main direction (x axis), and the 
vertical direction is the prong direction (y axis). 
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Fig. 2. Ground state wavefunction of a one-prong graph. We 
label the main direction x and the prong direction y. The com¬ 
ponent along the main direction has a kink due to the presence 
of the prong. 


damental limits without requiring any additional features such 
as donor-acceptor groups. Finite potentials will not lead to a 
phase disruption, which is quantified by a slope discontinuity 
of the wavefucntion at one point. The fact that even a small 
prong, which carries a small fraction of the electron charge, has 
a profound effect suggests that symmetry breaking in isolation 
cannot explain the observations. 

Fig. 1 shows a one-prong quantum graph which resembles 
a linear molecule with a side group normal to the main chain. 
The parts of the ground and excited state wavefunctions that lie 
along the main direction can be tuned to optimum shapes by po¬ 
sitioning a short prong to provide a pathway for charge to flow 
from the main charge-transfer axis, interrupting the flux along 
the main direction and creating a phase disruption in the wave- 
function along that direction. Fig. 2 shows the full ground state 
wavefunction xpQ (x, y) for the one-prong graph, where it is seen 
that the electron wavefunction along the main direction has a 
kink caused by the presence of the prong. We note that com¬ 
putation of quantum graphs requires solutions in both x and y 
directions [12] but the x direction is the dominant contributor to 
the response for this graph. 

Fig. 3 is the key result and shows how the behavior just de¬ 
scribed arises. The part (pnix) = ipn{x, 0) of the total wavefunc¬ 
tion xpn{x,y) lies along the main direction, and the first three 
levels (shown in thicker line type) are the most significant con- 



Position X along main direction 


Fig. 3. The amplitudes (pn(x) of the wavefunctions along the 
main direction are plotted for levels n = 0,1,2,3,4,5 for three 
prong lengths. The ground and first two excited states are plot¬ 
ted using thicker lines. The dashed line indicates the location 
of prong of length h, which is normal to the main direction. 

The components of the wavefunctions along the prong are not 
shown. The values of ^xxx and jxxxx are shown for each of the 
three prong lengths. 
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tributors to ^xxx for this graph topology. In all three panels, 
there is a component of the total wavefunction along the prong 
going into the page at the location of the dashed line. The length 
of the prong, h, is indicated on each panel. Also indicated are 
the values of ^xxx and jxxxx, again normalized to their maxi¬ 
mum values. 

The top panel of Fig 3. shows the effect of the phase disrup¬ 
tion caused by a very short prong. The wavefunctions are par¬ 
tially localized in such a way that the change in dipole moment, 
Xii is quite large when an electron is excited from the ground 
to the first excited state, but the overlap of the two wavefunc¬ 
tions is nearly zero, making xqi 0 and eliminating this exci¬ 
tation's contribution to the nonlinearity. Since this is the largest 
contributing term to ^xxx/ foe response is very small, despite 
the large change in dipole moment: Both large wavefunction 
overlap and a large dipole moment change are simultaneous 
requirements for a significant response. Similar remarks apply 
for jxxxx/ though the falloff with prong length is faster. Also 
shown are the next three wavefunctions to illustrate that they, 
too, have phase disruptions, but their contribution is small due 
to their much larger excitation energy. 

As the prong is lengthened, the middle panel shows that 
both a large change in dipole moment and good overlap be¬ 
tween the first two states persist, indicating large Xn and Xqi. 
Consequently, both the first and second hyperpolarizabilities 
rise significantly: A short, but not too short, prong results in 
ideal conditions for large response. Finally, as the prong is 
further lengthened, the shapes of the wavefunctions along the 
main direction approach those of a graph without a prong but 
with less charge along the main direction due to charge being 
drawn into the prong. Thus, the change in dipole moments 
upon excitation grows smaller, and ^xxx drops significantly. But 
Jxxxx is effected more so as it is more sensitive to the drop in 
changing dipole moments. 

We can easily identify the quantitative changes in the rele¬ 
vant transition moments as the prong length is tuned from a 
short value to a much longer one by examining the first few 
contributing terms to ^xxx from Eqn. (1). Fig. 4 displays the 
three lowest-energy contributing terms, their sum, and the full 
sum over all states for ^xxx as a function of the prong length 
h (top), along with the relevant transition and dipole moments 
as a function of h (bottom). For very short prongs, the change 
in dipole moment is largest for the first transition, but the 
wavefunction overlap Xqi these states is essentially zero. As 
the prong length increases, the change in dipole moment drops 
incrementally but remains large, while the wavefunction over¬ 
lap moves away from zero, thereby generating a significant con¬ 
tribution to the response. Note that for this graph, the three- 
level model for ^xxx is nearly exact. 

The analysis reveals the existence of an optimum prong 
length and position such that the overlap and the change 
in dipole moment between contributing states are both large. 
Short prong lengths of order 7% of the length along the main 
direction, located approximately 40% from either end, gener¬ 
ate a large first hyperpolarizability for the single electron graph, 
and the response drops dramatically as the prong length is de¬ 
creased to much shorter values but also as the prong exceeds 
this optimum value. 

These results generalize to a non-interacting, many electron 
system. Fig. 5 shows ^xxx as a function of the number of elec¬ 
trons Ne. The Fermi level is the highest occupied single electron 
state. The large variation in the response with electron count 
is a consequence of the specific geometry chosen for the graph. 



Prong length b 


Fig. 4. The full sum-over-states expression evaluated for ^xxx/ 
the three-level model sum + ^21 + ^12 + 

and contributions to ^xxx from the various three-level model 
terms are shown (top) with corresponding values of the rele¬ 
vant transition moments (bottom) as functions of the prong 
length. The dependence of the moments on prong length is 
as expected from the behavior of the wavefunctions functions 
in Fig. 3. (Reprinted with permission of Journal of Nonlinear 
Optical Physics and Materials.) 


For a given number of electrons, there exists a prong position that will 
maximize the response by optimizing phase disruption of the wave- 
functions for states near the Fermi level. For the graph studied, 
an electron count of Ne = 18 produces a large response. Fig. 6 
shows the wavefunctions along the main direction for the states 
directly below, at, and above the Fermi level with strong phase 
disruption for wavefunctions corresponding to levels nine and 
ten. Identical principles to those generating large responses 
in the one electron quantum graphs apply for many electrons 
when the electronic states nearest to the Fermi level show phase 
disruption, as is the case for the molecule in the inset of Fig. 5, 
which we propose as a prototypical large-system. 

The general principles delineated here reveal that electronic 
nonlinear optical molecular designs may be enhanced by judi¬ 
cious placement of a small side group along the main chain of a 
long molecule. Large nonlinearities necessarily arise when the 
most significant terms in the sum-over-states are enhanced by 
tailoring the shape of the ground and lowest excited state wave- 
functions through a phase disruption introduced by this side 
group when it is positioned closer to the center than to an end 
of the chain. Short side groups create significant phase disrup¬ 
tion while bleeding off very little charge from the main compo¬ 
nent of the dipole moment, and at the same time, maintain good 
wavefunction overlap for side groups that are not vanishingly 
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Fig. 5. ^xxx as a function of electron count Ne. The inset shows 
an 18-electron polyyne molecule whose geometry and topol¬ 
ogy matches the calculated one, giving ^xxx ~ 0-34. 

small. This behavior strongly supports the view that symmetry 
breaking alone is not the cause. These results also hold for the 
second hyperpolarizability suggesting that to any order in non¬ 
linearity promising molecules may be significantly improved 
by adding a side group off the main chain. 

One might argue that quasi-linear structures with infinite po¬ 
tentials at either end and a barrier or well potential in the cen¬ 
ter might give similar enhancement. Alternatively, a slant well 
modeling a donor-acceptor potential across a main chain can 
also lead to a large response[13] without the need of a slope dis¬ 
continuity in the wavefunctions. However, it is difficult to engi¬ 
neer a step, bump or slant well potential with a strong-enough 
potential energy gradient. In contrast, phase disruption relies 
only on adding a branch that changes the boundary conditions. 

In conclusion, we have shown that a strong phase disruption 
at a prong's attachment point significantly enhances the nonlin¬ 
ear optical response. Such structures are simpler to make than 
using end groups to tilt the potential. It is likely that molecules 
with ^xxx below the gap, designed using somewhat dated prin¬ 
ciples, do not exhibit the required wavefunction shapes to pro¬ 
duce large changes in the dipole matrix elements and, at the 
same time, maintain large spatial overlap between the con¬ 
tributing eigenstates. The quantum graph models that we re¬ 
port here, and the corresponding polyyne structure that we 
propose, should have large phase disruptions corresponding to 
states near the Fermi level, and thus an ultralarge hyperpolariz¬ 
ability. 

Funding. SMM and MGK thank the National Science Foun¬ 
dation (ECCS-1128076) for generously supporting this work. 
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